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Abstract 

The generalization of the Yang-Baxter equations (YBE) in the presence of Z2 grading 
along both chain and time directions is presented. The XXZ model with staggered 
disposition along a chain of both, the anisotropy ±A, as well as shifts of the spectral 
parameters are considered and the corresponding integrable model is constructed. The 
Hamiltonian of the model is computed in fermionic and spin formulations. It involves three 
neighbour site interactions and therefore can be considered as a zig-zag ladder model. The 
Algebraic Bethe Ansatz technique is applied and the eigenstates, along with eigenvalues 
of the transfer matrix of the model are found. The model has a free fermionic limit at 
A = and the integrable boundary terms are found in this case. 

This construction is quite general and can be applied to other known integrable models. 
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1 Introduction 



The construction of integrable chain models with inhomogeneous model parameters is very 
important in many aspects. Besides its own interest in mathematical physics, its possible 
applications in the problems of condensed matter physics is hard to miss-estimate. Even 
more, we believe that in order to proceed further in the non-critical string theory and 
cross the so called c = 1 barrier, it is necessary to consider strongly correlated fields 
on a line, instead of the free field degrees of freedom on the world sheet of the string (as 
in d < 1 case) . 

The crucial problem for strings in c? > 1 is the problem of the tachyon, which means 
a non-trivial structure of the ground state. Therefore, the investigation of strongly corre- 
lated chain systems, which can be considered on a fluctuating line, is actual. 

In this sense, in order to be able to take into account non-trivial ground state cor- 
rectly, it is straightforward to consider an integrable 2d model on a fluctuating surface. 
One attractive possibility is the consideration of inhomogeneously shifted spectral param- 
eters in i?-matrices along a chain. It is tempting to consider first the case of staggered 
inhomogeneity. 

An interesting model of hopping fermions in the Manhattan lattice, where the hopping 
parameters are staggered in both time and chain directions, appeared in connection with 
three dimensional Ising Model 0, In a particular case of hopping parameters, this 
produces a N x N numerical transfer matrix 0] in connection with edge excitations in 
the Hall effect. 

Since the pioneering work of H. Bethe and the essential input made by C. Yang 
the theory of 2(i-integrable models is well developed presently and has its formula- 
tion in technique, called Algebraic Bethe Ansatz (ABA) (or Quantum Inverse Scattering 
Method(QISM)) 0, ||. 

The basic element of this method is the matrix Rij{X, u) (A is the model parameter and 
u is the spectral parameter), the product of which along a chain defines the monodromy 
matrix T{X,u) as follows 

N 

T{X,u) = l[R^^{X,u). (1) 

i=i 

The trace of T{X,u) over states along the chain defines the transfer matrix t{X,u), the 
L's power of which (L is the lattice size in time direction) is the partition function of the 
model. 

When any two transfer matrices of the model with different spectral parameters com- 
mute then we have an infinite number of conservation currents, which causes the integrabil- 
ity. The sufficient condition of commutativity of the transfer matrices is the Yang-Baxter 
equation (YBE) for the i?-matrix, any solution of which defines an integrable model. 

Usually the integrable models under consideration are homogeneous along a chain, 
namely, the spectral parameters u in the product (|ID are equal. There were attempts to 
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construct inhomogeneous models, for example we can mention the successful attempt in 
the article of the model with alternating spin-i, spin-1 chain, but the main problem 
here is to obtain a local Hamiltonian (local conservation laws). It is obvious that, with 
the shift of the spectral parameter u in the product of Roj{X,u) matrices at the site j 
by arbitrary values Zj, we still fulfill the YBE (and hence the commutativity of transfer 
matrices) however the Hamiltonian of the model, and the other conservation quantities, 
will be local only if the shifts have periodic structure. Due to periodicity which goes far 
from neighbour sites, the local Hamiltonian will contain interactions between sites within 
the period and therefore can be represented as multi-leg ladder model. This type of ladder 
models differs from ones, constructed in [2^] by extension of the symmetry algebra, in 
2l| by construction first of the ground state, in [^] by use of higher order conservation 



laws and others P3l E4 



The aim of the present article is to construct an integrable model with commuting 
transfer matrices for different spectral parameters, where the i?-matrices in the product 
(H) are staggered by their construction as well as by shift of the spectral parameter u and 
where the Hamiltonian (and the other conservation laws) is local. The shift of the spectral 
parameter produces an additional model parameter 6. Similar model, where only the shift 
of spectral parameter was considered (no change of the structure of the i?-matrix), was 



constructed earlier in 10, 11 



We have generalized the YBE for that purpose in order to ensure the commutativity 
of the transfer matrices and demonstrate that they do have a solution in a case when the 
basic i?-matrix is the one of the integrable XX Z- mo del. 

In Section 2 the generalized YBE is formulated and the solution is found for the basic 
XXZ-model. 

In Section 3 we calculate the Hamiltonian for the XXZ model in both, fermionic 
and spin formulations, and demonstrate that it is local. Because of the invariance of 
the model under double space translations, the Hamiltonian contains interaction between 
three neighbouring points (instead of nearest-neighbour links in the ordinary XXZ). It 
can hence be represented as a zig-zag ladder Hamiltonian, where we have two terms, 
corresponding to SU{1, 1) Heisenberg models on each of two chains, and the interaction 
between chains terms of topological form, which are written on the triangles. In the free 
fermionic limit (A = 7r/2) we have obtained two independent sets of fermions, hopping 
only at the odd or even sites. 

This method, as one in [|l^, |l^], provides a way of construction of the integrable ladder 



models. 

Algebraic Bethe Ansatz (ABA) technique is applied to present model in Section 4 in 
order to find the eigenstates and eigenvalues of the transfer matrix. 

In Section 5 the open chain problem is considered and integrable reflection /C-matrices 
are found for the staggered XX-model. We also write the corresponding boundary terms 
in the Hamiltonian. 

We think that our approach for the construction of integrable models with staggered 
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model parameters is quite general and can be applied to other integrable homogeneous 
models. 



2 The definition of the YBE for the staggered model 
and its solution for the XX Z chain 

The key of integrability of the models is the YBE, which implies some restrictions on 
the Rij{X,u)-matnx, basic constituent of the monodromy matrix (|ID. The YBE ensures 
a local sufficient condition for the commutativity of the transfer-matrices t{u) = trT{u) 
at different values of the spectral parameter u, which corresponds to the rapidity of 
pseudoparticles of the model 

[Tiu),riv)]=0. (2) 

We use in the future the fermionization technique (alternative to Jordan- Wigner trans- 
formation), developed in [|12|, |T3|, Q (see also [jl5|) and work with the /^-operators (rather 



than matrices), expressed in Fermi-fields Cj, c^, (i = 1, ....iV is the chain site). 

By definition Raj acts as a intertwining operator on the space of direct product of the 
so called auxiliary Va (f ) and quantum Vj (u) spaces 

Raj iu,v) : Vaiu) ® Vjiv) ^ Vjiv) ® Vaiu) (3) 

and can be represented graphically as 



Raj{u,V) = 



I Vjiv) 
Figure 1: Raj matrix 



The spaces Va{u) and Vj{v) with spectral parameters u and v are irreducible represen- 
tations of the affine quantum group UgQ, which is the symmetry group of the integrable 
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model under consideration. Provided that the states | a) G K and | j) G Vj- form a basis 
for the spaces Va and Vj, following |12| we can represent the action of the operator Raj as 

Raj\ji)®\a^) = iRa,):f^\a2)^\j2), (4) 

where the summation is over the repeating indices 02 and j2 (but not over ai and ji). 
By introducing the Hubbard operators 

X:i=\a2){a^l X|^=|j2)(ji| (5) 
in the graded spaces Va and Vj correspondingly, one can rewrite (|[) as 

Raj = Raj I jl) I I (jl I = (RajTalf^ ' '^2) I J2) (cH | (jl | 

= i-lf^'^'^'^HRa^rafjlX^lXjl (6) 

where the sign factor takes into account the possible grading of the states | Oj) and | ji), 
p{ai) and p{ji) denote the corresponding parities and the summation over the repeating 
indices. 

In terms of operators Rij the matrix valued YBE can be written in the following 
operator form 

Rab{u, v)Raj{u, w)Rbj{v, w) = Rbj{v, w)Raj{u, w)Rab{u, v) . (7) 



We use for Vq and Vj the Fock space of the Fermi fields Q q, with basis vectors |0)i 
and for which 

' 1 - ni cf 

Ci Tli 



iX.)a 



we have a fermionization of the model, which is equivalent to Jordan- Wigner transforma- 
tion. 

In case of XX Z model, the fermionic expression of the operator Rij{u) is easy to 
obtain by use of formulas (||) and (||), and the standard expression for the R^a As a 
result one obtains 

Rjk{u) = a{u)[-njnk + {1 -nj){l - Uk)] +b{u)[nj{l - Uk) + {I - nj)nk] 

+ c{u) [c+Cfc + 4c,] . (9) 

Let us now consider Z2 graded quantum Vj-,p(f ) and auxiliary Va,a{u) spaces, p, a = 0, 1. 
In this case we will have 4x4 i?-matrices, which act on the direct product of the spaces 
Va^a{u) and V^-,p(f), (cr, p = 0,1), mapping them on the intertwined direct product of 
Va^a{u) and V^-,p(t,) with the complementary a = (1 — cr), p = (1 — p) indices 

Raj,ap (m, v) : Va,aiu) ® Vj,p{v) ^ Vj,p{v) ® K,<x(m). (10) 



'"^if the dimensions of the spases Vq and Vj are more than two, then one can use Fock spaces of more 
fermions, see [flSi 
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It is convenient to introduce two transmutation operations Li and 62 with the property 
'^i = '•i = '^'^ quantum and auxihary spaces correspondingly, and to mark the 

operators Raj,ap as follows 

Raj, 00 = Raj, Raj,01 = Raj, 

Raj AO = Raj, Raj,U = Raj^ ■ (H) 

The introduction of the Z2 grading in quantum space means, that we have now two 
monodromy matrices Tp,p = 0, 1, which act on the space Vp = Y[f=i ^j,p mapping it 
on Vp = nf=i Vj^-p 

Tp -.Vp^V-p, P = 0,1. (12) 

It is clear now, that the monodromy matrix of the model, which should define the 
partition function, is the product of two monodromy matrices 

T{\u) = To{\u)Ti{\u). (13) 

Now, because of the grading in the auxiliary space, we would like to construct the 
monodromy matrices Tq 1 as a staggered product of the Raj and R''^^ matrices. Let us 
define 

AT 

Ti(A,m) = \[Ra,2j-i{\u)Rl^^{\u) 

i=i 

N 

where the notation R in general means the different parameterization of the i?-matrix via 
models A and u parameters and can be considered as an operation over R with property 
R = R. 

Graphically the formulas ( [T^ ) can be expressed as in figure 0. 

In order to have a integrable model with commuting transfer matrices (|I3D for different 
spectral parameters 

[trT(A,u), trT{X,v)] = (15) 
it is enough to have the following relations for the Tcr{u) = trTcr{X,u), {a = 0, 1) 

T„{\,u)ti^„{\,v) = f^{X,v)fi^a{\u), (7 = 0,1 (16) 

It is not hard to see, that in order to ensure the commutativity condition (p!5| ) the 
R- and ^-matrices in the expression (|14|) should fulfill the following two Yang-Baxter 
equations 

Rabiu,v)R'^]^iX,u)Rk,,iX,v) = Rl]^iX,v)Ra,jiX,u)Ratiu,v) (17) 

5 




Figure 2: Monodromy matrices Tq and Ti 

and 

Rat{u,v)K%,{\,u)Rl-^^^,iX,v) = Rl%,{\v)R'l.^,{\u)K,{u,v). (18) 

By inserting the expression (P) for the i?-operator into (p!7|) one can obtain 20 equations 
on parameters a, 6, c and a, 6, c, which after defining the transposition operation ii and 

~ as 

a'^{u) = a{u), b'^{u) = -b{u), c'^ = c{u) 

d{u,v) = a{u,v), b{u,v) = —b{u,v), c{u,v) = c{u,v), (19) 

are reducing to following three equations 

a{u, v)c{u)a{v) = c{u, v)a{u)c{v) + b{u, v)c{u)b{v), 
b{u, v)d{u)c{v) + c{u, v)c{u)b{v) = a{u, v)b{u)c{v), 
c{u, v)b{u)a{v) = b{u, v)c{u)c{v) + c{u, v)d{u)b{y). (20) 

It is easy to recognize in equations (|20| ) the ordinary YBE for the XX Z model, where 
we suppose that the parametrization of the variables d{u), b{u), c{u) differs from a{u), 
b{u), c{u). We are going to introduce a new model parameter here for them. 
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The second set of YBEs (|T8|) gives 

a{u,v)c''^{u)a'-^{v) = c{u,v)a''^{u)c''^{v) - b{u,v)c''^{u)b'-^{v), 
-h{u,v)a'''^{u)c'''^{v) + c(m, t;)c''^(u)6'^(f) = a{u,v)h'''^{u)c''^{v), 
c{u,v)h'\u)a'\v) = -b{u,v)c''{u)c'"'{v) + c{u,v)a'%u)V^{v). (21) 



which differs from the first set (pOD by the transposition i2 of the variables, by displaced 
bars from the a{u), b{u), c{u) on to a{v), b{v), c{v) and change of sign b{u,v) — > —b{u,v). 
This last change of the sign comes from the fact that, in the first set of YBEs ([T7|), the 
variables b{u, v) stationed in the right side of equations, whilst in the second set (|T8|) - in 
the left side. 

The solution of the equations ( pOD is the same as for ordinary XX but we will make 
a different shifts of the spectral parameters for the a{u), b{u), c{u) and a{u), b{u), c{u). 
Namely, if we take 

a{u) = sin(A + u), a{u) = a{u) = sm{X + 9 — u) , 

b{u) = sin(M), b{u) = b{u) = sm{9 — u), 

c{u) = sin(A), c{u) = c{u) = sin(A), (22) 

then the solution of (pOD for a{u,v), b{u,v), c{u,v) is 

a{u,v) = a{u — v) = sm{X + f] — u — v), 
b{u,v) = b{u — v) = sm{f] — u — v) , 

c{u,v) = c(u — f ) = sin(A), (23) 

where u = 9 — u and 9 is the new parameter of the model. 

It is easy to see after simple calculations, that the second set of YBEs ( |21| ) is compatible 
with the first set (^OD by producing the same solution (pSf) , if we define the i2 transposition 
as follows 



a'^(M) = sin(A-M), {u) = b{u) , c'^ = c{u), 

a'^{u) = sin(A-M), b'^{u)=b{u), c'' = c{u). (24) 

Therefore, the formula (|22|) together with the definitions (0) and (p^ defines a so- 
lution of the double set of YBEs (0) and ([T8|). The new parameter 9 appeared in the 
model in addition to A (or A) in the ordinary XXZ-model. 

As it is well known, the variable 

A = —-^ , , — ^ = cos A = const 25 

characterize the anisotropy of the product of the spins (Tj and CTj+i in z-direction in the 
Heisenberg model. 
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It is important to mention here, that the transpositions Li and L2 changes the sign of 
A, = At2 = —A. This means, that in a definition (|I^) of the monodromy matrices 
To,i(A,u) we have a product of i?-matrices with alternating anisotropy parameters ±A. 



This fact essentially differs our construction from the one in the articles |T^, |rT|, where 
the simple shift of the spectral parameter is considered. 



3 The Hamiltonian of the model. Ladder form 

The prove of existence and the calculation of the local Hamiltonian in our model is easy 
to carry out technically in a so called check-formalism, where instead of usual i?-matrix, 
we consider the matrix R = PR, [P being the permutation map P = x®y^y®x). In 
matrix elements 

(26) 



.a'b' 
ab 



R. 



b'a' 
ab 



Instead of the operator Raj in (^) let us define 



Ri 



(R^^r^b' \a')\b'){b\{a\={-l) 



p{a')ip(b')+p{b)) 



(R 



\a'b' -ya -yb 
ijJab ^a'^b'^ 



(27) 



which after the fermionization (M) will have the form 



Rjk{u) = a{u) [rijUk + (1 - nj){l - Uk)] + c{u) [nj{l - Uk) + (1 - nj)nk] 
+ h{u) [cjck + 4cj] . 



{2i 



and where the parameters a(u), h{u), c{u) are defined by formulas (^21) . 

The monodromy matrix (|n[) in this formalism has the following expression 

N 

,2i+l (A, u)R2j+i^2j+2{\ 

N 

To{X,u) = ^^i+l,2,+2(A,^)i?2^/|2,2,+3(A,«)■ 
i=l 



(29) 



It is convenient to represent graphically the fermionic ^,j+i-operators as in figure 0. 




Figure 3: Ri 



8 




Figure 4: Monodromy matrix 



Then, the full monodromy matrix T(A, u) = Tq{X, u)Ti{X, u), defined by use of formu- 
las (^), can be pictured as in figure ^, which also demonstrates the order of writing of 
the operators in the product. 

In order to be able to write a local Hamiltonian, it is necessary to have a value 

Rij{uo) = lij. (30) 

Then the logarithmic derivative of the transfer matrix at this point will define a local 
Hamiltonian 

u _ d\ntrT{u) , 

-n — \u=uo ■ l-Jj-j 

Technically, in order to find a Hamiltonian, one should expand -Rj^-operators around 
Mo and extract the first order terms over (u — uq) in the product T(A, u) = To(A, m)Ti(A, u). 
In ordinary homogeneous case by inputting the first order expansions of the i?-operators 
around unity into the product 

AT 

T(A,M) = n4.+i(A,w) (32) 
i=i 

it is easy to see that there is no scattering between long distance fermions and only the 
nearest-neighbour hopping terms are contributing into the Hamiltonian. 

In the present model we have a two different -R-operators in the product (P^), only 
one of them becoming unity for some value of the spectral parameter. Let us take Mq = 0. 
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Then, by expanding around the point u = the R, R'''^'''^, R , R operators, defined by 
formulas (|T|) , , and (^), obtains 

R2i,2i+l{u) = sin Xl2i,2i+1 + m(cOS XP2i,2i+l + C2iC2i+l + C^i+iC2i), 

-R2i'2i+l(«) = sin \l2i,2i+l - m(cOS \P2i,2i+l + cfci+i + cf^iQ) , 

^2i-l,2i(^) = ^2i-l,2i(^) - U [C0S(A + e)P2i-l^2i " COS 6 {c^^_^C2i + 4iC2i_l)] , 

^2Ll,2i(«) = R2i-l,2i{.-0) + M [C0S(A - 9)P2i-l,2i - COS r7(c^._iC2i + C^jCsi-l)] , (33) 

where = riirii^i + (1 — ni){l — rii^i) and R2i-i,2i{6) defined by formula (^) with 

a(e) = sin(A + ^), 6(0) = sin ^, c(0)=sinA. (34) 

It is easy to see in the Figure | that, although the ^2i-i,2j-operators are not unity 
at u = 0, the Hamiltonian will still be local. The scattering due to ^2i-i,2roperators 
will induce an interaction involving at most 4 fermions. For example, the excitation 
Hamiltonian 7^23 for the fermions 2 and 3 (see Figure ^ due to scatterings in RuiO) and 
R34,{9) may create the interaction between 1, 2, 3, 4 states only. 

Now, by inserting the expressions (p3D and ( ^If ) for R , R , R and R''^''^ into the 
monodromy matrix T{u) = Tq{u)Ti{u) , in the linear approximation of u we obtain 

T{u) ^ (sin^ A sin(A + 9) sin(A - 9))^ + ^Yl ' (^^^ 

with 

sin A sin(A + 9) sin(A - 9)Hj = 

+ (-!)■' sin6'sin2A(l - rij^i - ?2j+2)(c+Cj+i - c+^_iCj) 

+ (-1)-'^^ sin^ 0COS A [njnj+2 + (1 - ^i)(l - "-^+2)] 

+ sm9 [(sin(A + {-iy9){l - uj+i) + sin(A - (-l)^'0)n,+i) c+c,+2 

- (sin(A - {-iy9){l - Uj+i) + sin(A + {-iy9)n,+i) c+^^c,] . (36) 

In the limit A = 7r/2, which corresponds to XX-model, we will have two types of 
independent free fermions, hopping separately at the even and odd sites of the chain with 
the Hamiltonian |^ 

Hj = tan^ 9 [c+9+2 - 4_^2Cj] ■ (37) 

It is convenient now to write the Hamiltonian (|36D in a ladder form represented graph- 
ically as in Figure ^. 

^It seems to us that the model considered in do not reproduces the free fermionic hmit at A = 7r/2, 
as in case presented here. 
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2 J + 1 2j + 3 




Figure 5: Zig-zag ladder chain 

Let us consider the even (2j) and the odd (2j + 1) points of the chain as a sites {j) 
of two different chains labeled by s = and 1 correspondingly. The Fermi fields will be 
marked now as 

Cj^s = C2j+s, s = 0, 1 (38) 

With these notations it is straightforward to obtain from the expression dH) the 
following ladder Hamiltonian 

sin(A + 6) sin(A - 6) sin(A)ifj- ^ = 

(-l)''sin6'sin2A (c+^Cj-^+i - c+^+iCj-^) {uj^i^s+i - 

+ (-1)''+^ sin^ 9 cos X[nj^snj+i,s + (1 - nj^s)^^ - nj+i^s)] 

+ sine (sin(A + (-1)'^)(1 - + sin(A - ,+i) 

-4+i,sCj,s (sin(A + (-l)^^)n,,,+i + sin(A - {-iye){l - n,-,+i))] . (39) 

Hence, we have obtained a simple integrable ladder Hamiltonian, which has a free fermionic 
limit. 

We can also write down the Hamiltonian of the model in a spin language, calculated 
by use of the ordinary representation for i?"^,* as 

/ a \ 

pa'6' ^ C 6 

6 c ■ 
\ a / 
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In ladder form the result looks as follows 

sm(A + 6) sin(A - 6) sm{X)Hj = 
(-lY 

sin^ 9 cos A [^^^l+i,^ + " ^s^^l+i, J 

+ sin e sin A {cos A(-l) (a+.+ia^+i,, - cr J s+i(^^+i,s) 
{-ly cos A {(r^,s(^ls+i - ^Is^t+i) 4+1,^ 

cos 9 - (^J,s^ls+i^J+i,s} } (40) 

As one can see from this expression, the first term gives XXZ models for the each 
leg of the ladder with anisotropy parameters A = — 1 and therefore can be considered as 
SU{1,1) Heisenberg chains. The other terms, which represents the interaction between 
chains, are written for the each triangle of the ladder (see Figure ^ and has a topological 
form 

e crj-,sCrj-,s+iO-j+i,^, (41) 
where e"'"^ is an anisotropic antisymmetric tensor, defined by the coefficients in the formula 



4 Algebraic Bethe Ansatz (ABA) solution of the model 

The technique of ABA , called also Quantum Inverse Scattering Method (QISM), 
essentially was developed in the works of Baxter ^ and Faddeevs group |0 . 

In order to carry out ABA it is convenient to work in conventional (not braid) for- 
malism and use the YBE for i?-operators in the form of formulas ( p!7| , |T8|). Let us first 
define the L-matrix as 

{L,{u)):, = {a I R,, I a% = {-lyi'^'M"') {R,,{u)t,, 

_ / a{u){l — rij) +h{u)nj c{u)c^ \ 
\ c{u)cj —a{u)nj + b{u){l — rij) J ' 

which is a matrix in the horizontal auxiliary space with operator value entities and acting 
on quantum space Vj. The matrix elements between auxiliary states | a') and {a \ of the 
monodromy operators To,i(A,u) defined in ( p!^ 

can be expressed as a product of the L-matrices as follows 

iTi):i = {L{)Z{u){v,r:M {l\xz-^{u), 

(To)Z = {L[TAn){Lr)Z{u) (Lr):--(«). (42) 
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By use of equations ( p!7| , p!8| ) one can obtain the graded YBE for the monodromy 
matrices T^, s = 0, 1 as follows 

(_i)(P(.')-i-p(."))PK)^.?^,(,,,) (To:: («) (To):: (.) 

= (-l)^(^')(^('^')+^(^)) {T,t, (v) (To):, iu)R^Uu,v). (43) 
Now, following the procedure of ABA, let us define empty fermionic state 

2N 

I ^) = n 1 (44) 

i=l 

as a test vacuum of the model and demonstrate that it is the eigenstate of the transfer 
matrix t{X,u) = ti{\, u)tq{X, u) = strT{X,u). 

The expression (|4^) for L-matrix shows that his action on | 0)j produces upper trian- 
gular matrix, therefore the action of the monodromy matrix (Tg)'^ on | defined by the 



formulas (^2]) and (^^ , will also have upper triangular form and can be calculated easily 

as 

Now it is obvious, that | Q) is eigenstate of r(A, u) with eigenvalue 



u,{u) = [a{u)a''{u)f -[b{u)b'''{u)f. (45) 
One can see from the expressions (H3|) that the operators Cs{u) act on | Q) as the annihi- 



lation operators, while Bs{u) act as the creation operators. That is why it is meaningful 
to look for states 

\vi,V2,...Vn)o = Bi{vi)Bo{v2)...B^{vn) \^)x, X = n{mod2) , (46) 

as an n-particle eigenstates of r(n) with spectral parameters Vi, i = 1, ...n. In order to 
check whether this is true we do not need to have an exact form of operators Bs{u), we 
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need only to know the algebra of operators As{u), Ds{u) and Bs{u), which can be found 
from the YBE (]43| ) as follows 



and 



Ai{u)Bo{v) 

Do{u)Bi{v) 
Di{u)Bo{v) 



a[v 


-u) 


h{y 


-u) 


a{y 


~u) 


' h{v 


~u) 


a{u 


-V) 


h{u 


-v) 


a{u 


-v) 



h{u — v) 



Bi{v)A^{u) - 



Bo{v)Di{u) + 
B,{v)Do{u) - 



b{v 


-u) 


c{v 


-u) 


b{v 


-u) 


c{u 


-v) 


b{u 


-v) 


c{u 


-v) 


h{u 


-v) 



■Bi{u)A^{v) 

B^{u)Di{v) 
B,{u)Doiv) 



(47) 
(48) 

(49) 
(50) 



The first terms in the right hand side of equations (^^ and (^9[) are so called "wanted" 
terms and they are producing the eigenvalues z/(m, fi,f2, ..f„) of the state | fi,f2, ■■■Vn)o 
as 

(-l)t 



I/o(M,fl, ■■■,Vn) 



l^l{u,Vi, ■■■,Vn) -- 



n n 

Y[ a{vi - u) [a'^ {u)a'''^ {u)f - JJ a(M - v,) [b'' {u)b'''^ (u)]^ 

1=1 

nr=i - u) 

n n 

Yl a{vi - u) [a{u)a"' {u)f - JJ a(M - Vi) [b{u)b'' {u)f 



(51) 



i=l 



But in order the state | fi,f2, ■■■Vn) to be an eigenstate of r{u) we need the cancellation 
of so called "unwanted" terms, produced by the second terms in the right hand side 
of equations (|47|) and (^91). This gives us the restrictions on the spectral parameters 
Vi,V2, ■■■Vn in a form of Bethe Equations (BE) 



a{Vj)a^^{v 



N 



n 



a[Vj - Vi) 



a[v^-Vj, 



j = l,...n. 



(52) 



These equations are similar, but differ by some sign factors from the Bethe equa- 
tions obtained for the model defined in |]T0| caused by the alternating anisotropy model 
parameter ±A along the chain. 
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The calculation of the eigenvalues of the Hamiltonian ( |39[ ) is straightforward by taking 
of logarithmic derivatives of (|5T|) at m = , which gives us 

" / 1 1 \ 

\sint>j sin(A + f j) sm{vi — 6) sm{X + Vi — 6) ) 

sm29 

^ sin(A + 6) sin(A - 9) ' ^ ' 

5 Open chain Hamiltonian and boundary terms of 
the model in free fermionic case (A = 7r/2) 



In this section we follow the approach and the equations of [|lg, |T^, [T8|, |T9[ to con- 
struct the reflection /C-matrix for our model in free fermionic case and calculate integrable 
boundary terms for the Hamiltonian. 

In ordinary integrable models one defines the double row transfer matrix as 

r{u) = tio {K+{u)T{u)K^{u)T{-u)-') 

X K^{u)Ri2{u)R2z{u)---Rl-i,l{'^)Rlo{u)) , (54) 
where N is the number of chain sites and C,{u) is defined by the unitarity condition 

R{u)R{-u) = C{u) = cos^(u) (55) 

for the XX-model. 

The reflection matrices K{u) (for the right side) and K~^{u) (for the left side) have to 
obey the reflection equations, the analogues of YBE on the boundaries, in order to ensure 
the commutativity of the double row transfer matrix ( p^ for different spectral parameters. 
This, together with the YBE equation, is a sufficient condition for integrability. The 
reflection equations read 

Ruiu - v)K^{u)Ru{u + v)K^{v) = K^{v)Ri2{u + v)K^{u)R2i{u - v) (56) 

and 

Rui-u + v)K+{u)Ri2{-u -v + n)Kl{v) = 

= K+{v)Ri2{-u -v + 7r)Kt{u)Ru{-u + v) . (57) 

After the fermionization (|D,(|D,(H), the solutions of the equations (p6D and (|57| ) have 
the following form 

K{u) = sin(^ + u)n + sin(^ — «)(! — n) (58) 
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and 

K^{u) = - sm(^+ - u)n + sin(^+ + n), (59) 

where ^ and are arbitrary parameters. 

In the present model, since we aheady have a two-row transfer matrix in the closed 
chain case, the open chain transfer matrix contains four rows, with two- normal and two- 
backward directions. 

Let us define 

t{u) = To{u)n{u) = C W^tro [K+{e - u)Toie - u)Koie - u)To)-\-e + u)] x 

X tn[Ktiu)Tl'{u)K,{u)T^^)-\-u)], (60) 

where -^'0(^^)5 -^1(^)5 Kq{u) and K^{u) are defined by equations (^) , (^9]) , while To(m) 
and Ti{u) are defined by (^) after the shift of the spectral parameter. 

Let us emphasize here, that because the boundaries of the Tq{u) and Ti{u) in (p9|) 
can be defined in a various ways, namely, the beginning and the end of the chain can be 
translated by one i?-matrix, we will have a variety of boundary transfer matrices. We 
now consider only one possible solution. 

The boundary transfer matrix t{u) can be expressed graphically as in Figure |^. 




Toi9-u) 
Figure 6: Boundary transfer matrix 



This way of writing the boundary transfer matrix t{u) is not convenient for the calcu- 
lation of the Hamiltonian, particularly for the bulk part, but it is convenient to analyze its 
commutation properties for different spectral parameters. It is clear that the boundary 
YBEs (reflection equations) are the same as (|56| ) and (^Tf ) for the ordinary case. The 
expressions (^H|) and (pD|), after the shift of the spectral parameter u ^ 9 — u for the 
row will ensure the commutativity of the boundary transfer matrices (^). Moreover, by 
use of YBE ( p!7| ) and ([TB|) we can translate the Tq"^(— 6* + u) row in Figure ^ to very top 
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by commuting with T(^(m) and (T^^) ^{—u) rows and bring it to the form, convenient for 
the calculation of the Hamiltonian: 



t{u) = tr [lC+{u)To{e - u)Ti{u)lC{u)T{\-u)T^\-e 



u 



(61) 



Here the effective reflection matrices }C^{u),IC{u) have been introduced and are defined 
as follows 



/C(n) = R{2u + e)Ki{u)R{e)Ko{e -u), 
/C+(n) = -R{e -2u)K^{u)R{-e)K^{e -u) 



(62) 



It is clear, that the boundary transfer matrix r(n) can be represented graphically as 
in Figure 0. 



t(u) 



K+(u) 





Ki(u) 




To{e-u) To{e-u) 

Figure 7: Identity relation for boundary transfer matrix 

The equation in Figure |^ is a consequence of the reflection equations (|56| ) and 
Now again, in order to calculate the boundary terms, we should input into the ex- 
pression of the boundary transfer matrix (|6lD, as in the case of the bulk Hamiltonian, 



the linear expansions of the operators R, R, defined by the formulas (^ - |3J) and the 
following linear expansions of the reflection i^-matrices around u = 



Ki{u) = sin^ + u cos ^[ui - {1 - Hi)] , 
Ki{e-u) = [sm{^ + e)ni + sm{^-e){l-ni)] + 

and extract the terms linear in v. After some calculations one can obtain 

N 



(63) 



(64) 
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where Hj is the bulk Hamihonian, which have been calculated in the previous section and 
given by the formula (^), while the right -Hr and the left- boundary terms have the 
following expressions 



H 



R 



+ 



+ 



tan^ 
1 



sin(e - e) 

.sin(e + e) 

sin(^ - 9) 
cos 9 \ [sin(^ + 9) 
cot C, 
cos^ 9 



+ 1 



2 tan 6* 



cot(^ + 



rtN + - 



sin(^ + I 
sin(^ — I 

cot ^ 
cos^ 9 



+ 1 



Cj^_2Cn 



sin(^ + 



sin(^ 



1 + 2tan6' 



- cot(^ 



(1 - Un) 



Hr 



— tan^ 9 cot ^n7v-i + tan^ 9 cot ,^(1 — r;,Ar_i) — 2 tan 6', 
sin 9 



sin ^ cos 29 



const 



(65) 



It is important to mention that, besides the translation of the boundary of the chain, 
there is another possibility for a modification of the boundary transfer matrix (^). We 
can take the return paths Tq'^{—9 + u) and {T[^)^^{—u) not above the direct chain, as it 
was considered and presented in the Figure ^, but below. Then it is not hard to see that 
this will interchange the left and right boundary terms (|65|) in the Hamiltonian. Finally, 
there are four integrable models with different boundary terms, only one of which being 
explicitly written here. 
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